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ABSTRACT 



We present the results of a new study of the molecular gas associated with the 
"extended green object" (EGO) G35.03+0.35. This object, very likely a massive young 
r"| , stellar object, is embedded in a molecular cloud at the border of an HII region. The 

Q^ 1 observations were performed with the Atacama Submillimcter Telescope Experiment 

(ASTE) in the 12 CO and 13 CO J=3-2, HCO+ J=4-3, and CS 3=7-6 lines with an 
angular resolution about 22". From the 12 CO J=3-2 line we discovered outflowing 
^ 1 activity of the massive young stellar object. We obtained a total mass and kinetic 

energy for the outflows of 30 M and 3000 M Q [km s^ 1 ] 2 (6 x 10 46 ergs), respectively. 
We discovered a HCO+ and CS clump towards the EGO G35.03+0.35. The detection 
of these molecular species supports the presence of molecular outflows and a dense 
molecular envelope with temperatures and densities above 40 K and 6 x 10 6 cm~ 3 , 

■ respectively. Using public near- and mid-IR, and sub-mm data we investigated the 
I spectral energy distribution confirming that EGO G35.03+0.35 is a massive young 

QO ■ stellar object at the earliest evolutionary stage (i.e. a class I young stellar object). 

■ By anlysing radio continuum archival data we found three radio sources towards the 

object, suggesting the presence of several young stellar objects in the region. Our radio 
continuum analysis is consistent with the presence of at least one ultracompact HII 
region and an hypercompact HII region or a constant-velocity ionized wind source. 
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C« t 1 INTRODUCTION proved that the 4.5 /im emission is due primarily to lines of 

molecular hydrogen, which are collisionally excited. 

Cyganowski et al. (2008) identified more than 300 Galactic 

extended 4.5 /im sources, naming extended green objects EGO G35.03+0.35 (hereafter EGOg35) is embedded in 

(EGOs) or "green fuzzies" , for the common coding of the [4.5 a dense molecular cloud at the distance of 3.5 kpc (Petriella 

fim] band as green in three-color composite Infrared Array et al. 2010) located on the border of the infrared dust bubble 

Camera images from the Spitzer Telescope. According to the N65 (Churchwell et al. 2006, 2007). According to Petriella 

authors, an EGO is a probable massive young stellar object et al. (2010) there are several young stellar object (YSO) 

(MYSO) driving outflows. The extended emission in the 4.5 candidates around N65, being EGOg35 the most prominent 

/im band is supposed to be due to H2 {v = — 0, S(9,10,ll)) source. Using IR and sub-mm fluxes measured from this 

lines and CO (y = 1 — 0) band heads, that are excited by the EGO, the authors performed an spectral energy distribu- 

shock of the outflows propagating in the interstellar medium tion (SED), showing that this source is indeed a massive 

(see Noriega-Crespo et al. 2004; Marston et al. 2004; Smith stellar object at the earlier stages of evolution with outflow- 

& Rosen 2005). Recently, De Buizer & Vacca (2010) reported ing activity. The 4.5 /im emission of this EGO has a bipolar 

the first direct spectroscopic identification of the origin of morphology, with one lobe to the NE and the other to the 

the 4.5 /im emission towards two EGOs using NIRI on the SW. The source presents maser emission of several molec- 

Gemini North telescope. In one of the observed EGOs, they ular lines (see e.g. Forster & Caswell 1989; Caswell et al. 

1995; Kurtz & Hofner 2005). Recently, CH3OH maser emis- 
sion at 6.7 and 44 GHz has been detected and analysed by 

* E-mail: sparon@iafe.uba.ar Cyganowski et al. (2009). The 6.7 GHz CH3OH masers are 
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Figure 1. Left: Spitzer-IRAC three color image (3.6 /jm = blue, 4.5 (im = green and 8 fim = red). The dashed rectangles are the regions 
mapped in 12 C0 J=3— 2 and HCO"*" J=4-3 (the largest one) and in 13 CO J=3-2 and CS J=7— 6 (the smallest one). Right: zoom in of 
the surveyed region. 



radiatively pumped by IR emission from the warm dust as- 
sociated exclusively with massive YSOs, while the methanol 
maser at 44 GHz is collisionally excited in molecular out- 
flows, and particularly at interfaces between outflows and 
the surrounding ambient cloud (see Cyganowski et al. 2009 
and references therein). 

Considering that EGOg35 is indeed a massive YSO 
evolving in a dense molecular medium, it is important to un- 
derstand the physical characteristics of its ambient medium. 
We investigated this ambient through several molecular lines 
obtained with the Atacama Submillimeter Telescope Experi- 
ment (ASTE). The results of this investigation are presented 
in this work. 



2 OBSERVATIONS 

The molecular observations were performed on July 14 
and 15, 2010 with the 10 m Atacama Submillimeter Tele- 
scope Experiment (ASTE; Ezawa et al. 2004). We used the 
CATS345 GHz band receiver, which is a two-single band 
SIS receiver remotely tunable in the LO frequency range of 
324-372 GHz. We simultaneously observed 12 CO J=3-2 at 
345.796 GHz and HCO+ J=4-3 at 356.734 GHz, mapping 
a region of 70" x 80" centered at RA = 18 h 54 m 0.7 s , dec. 
= +02°01'18.9", J2000). We also observed 13 CO J=3-2 at 
330.588 GHz and CS J=7-6 at 342.883 GHz towards the 
same center mapping a region of 40" x 50". The mapping 
grid spacing was 10" and the integration time was 60 sec per 
pointing in both cases. All the observations were performed 
in position switching mode. We verified that the off position 
(RA = 18 h 53 m 55.8 s , dec. = +02°07'49.5", J2000) was to be 
free of emission. 

We used the XF digital spectrometer with a bandwidth 
and spectral resolution set to 128 MHz and 125 kHz, respec- 
tively. The velocity resolution was 0.11 km s" 1 and the half- 
power beamwidth (HPBW) was 22", for all observed molec- 
ular lines. The system temperature varied from T sys = 150 
to 200 K. The main beam efficiency was 77 m b ~ 0.65. The 
spectra were Hanning smoothed to improve the signal-to- 
noise ratio and only linear or/and some third order poly- 
nomia were used for baseline fitting. The data were reduced 
with NEWSTAR and the spectra processed using the XSpec 
software package. 



To complement the molecular data we use near- and 
mid-IR and radio continuum data from public databases 
and catalogues, which are described in the corresponding 
sections. 



3 RESULTS AND DISCUSSION 

The source EGO G35.03+0.35 (EGOg35) is located at the 
border of an HII region which is delineated mainly by the 8 
^m emission. Figure 1 (left) shows a composite three-color 
image of a region towards EGOg35. The image displays three 
Spitzer-IRAC bands: 3.6 fim (in blue), 4.5 fim (in green) 
and 8 /im (in red). EGOg35 is the green structure inside the 
dashed rectangles, which represent the 70" x 80" and 40" x 
50" regions mapped in the molecular lines as described in 
the previous section. A zoom in of the surveyed region is 
shown in Fig. 1 (right). 

3.1 Molecular lines 

To study the molecular ambient where EGOg35 is embed- 
ded, we analyse the 12 CO and 13 CO J=3-2, HCO+ J=4-3 
and CS J=7-6 transitions, tracers of outflows and dense gas. 
Figures 2 and 3 display the molecular lines spectra obtained 
towards EGOg35. Most of the spectra are far of having a 
simple Gaussian shape, presenting asymmetries, probable 
absorption dips, and spectral wings or shoulders, which sug- 
gest that the molecular gas is affected by the dynamics of 
EGOg35. In what follows we study the gas kinematics. 

Figure 4 shows the spectra obtained towards the central 
position of EGOg35, that is the (0,0) offset in Figs. 2 and 3. 
In the case of the HCO + J=4-3 and CS J=7-6 lines, their 
brightness temperatures were scaled with a factor of x3. 
Figure 5 displays a zoom- in of the central 12 CO and 13 CO 
J=3-2 profiles in the intensity range that goes from —0.35 
to 2 K, in order to note some weak features. In the case of 
the 12 CO line can be appreciated a weak component (~ 5 
times above the rms noise) centered at ~ 41 km s . The 
parameters determined from Gaussian fitting of these lines 
are presented in Table 1. T m f, represents the peak bright- 
ness temperature, Vlsr the central velocity referred to the 
Local Standard of Rest. Errors are formal la value for the 
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Figure 2. 12 CO J=3-2 and HCO+ J=4-3 spectra obtained to- 
wards EGOg35. 



model of the Gaussian line shape. All the lines were well fit 
with more than one Gaussian function, which very likely in- 
dicates the presence of several components or/and spectral 
wings, usually signatures of outflows. In Table 2 we present 
Vmax and V m in, indicating the respective total widths of 
the spectra considering all the components, and J T m bdv, 
the intensity integrated over the whole profile. 

The 12 CO 3=3-2 spectrum presented in Fig. 4 shows a 
double peak structure like with a main component centered 
at ~54.2 km s _1 and a less intense component centered at 
~48.7 km s _1 . The 13 CO J=3-2 spectrum presents similar 
features. We conclude that these lines present an absorption 
dip at ~51.5 km s _1 which separates both mentioned veloc- 
ity components, showing that the lines are self-absorbed as 
it is usually found towards star-forming regions (Johnstone 
ct al. 2003; Buckle et al. 2010; Ortega et al. 2010). The ve- 
locity of the mentioned 12 CO and 13 CO dip is coincident 
(within the errors) with the central velocities of the SiO (5- 
4), H 13 CO+ (3-2), and CH 3 OH (5 2 ,3-4i, 3 ) lines observed 
towards the center of EGOg35 by Cyganowski et al. (2009) 
and with the CS 3=7-6 main component reported in this 
work. Therefore, we conclude that v ~51.5 km s _1 is the 
velocity of the ambient gas, and the other reported compo- 
nents (see Table 1) may be related to outflows or/and high 
velocity material. In Section 3.2 we analyse this contention. 

Regarding the HCO + J=4-3 emission, the profiles to- 
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Figure 3. 13 CO J=3-2 and CS 3=7-6 spectra obtained towards 
EGOg35. 
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Figure 4. Spectra towards the center of the observed region. The 
HCO + J=4— 3 and CS J=7-6 lines were scaled with a factor of 
x3. The vertical axis is main beam brightness temperature. 



wards the center of the surveyed region present a dip at v 
~ 54 km s _1 (see Fig. 2 bottom, and Fig. 4). Cyganowski 
ct al. (2009) presented a single point observation of H 13 CO + 
J=4-3 towards EGOg35. They report that the H 13 CO+ 
J=4-3 emission peaks at ~ 53.1 km s _1 with a FWHM 
Av ~ 5.4 km s" 1 . Taking into account that the H 13 CO + 
emission is optical thinner than the HCO + emission and 
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Figure 5. Zoom-in of the 12 CO and 13 CO J=3-2 profiles pre- 
sented in Fig. 4 in the intensity range going from —0.35 to 2 K. 



Table 1. Observed parameters of the molecular lines shown in 
Figure 4. 



Emission 








(K) 


(km s-i) 


12 CO J=3-2 


0.80 ±0.20 


41.12 ±0.56 




4.55 ±0.25 


48.68 ±1.25 




12.40 ±0.75 


54.20 ±1.20 




4.42 ±1.50 


56.90 ±1.20 




3.65 ±0.85 


59.10 ±1.25 


13 CO J=3-2 


4.40 ±0.45 


50.10 ±1.60 




12.85 ±1.70 


53.50 ±1.40 




2.48 ±0.60 


56.75 ±1.50 


HCO+ J=4-3 


3.10 ±0.8 


52.55 ±0.75 




2.60 ±0.8 


56.05 ±0.85 


CS J=7-6 


2.85 ±0.35 


53.15 ±0.65 




1.45 ±0.60 


57.05 ±0.50 


Table 2. Derived parameters 


of the molecular lines sh 


Figure 4. 






Emission 


~V min 






(km s-i) 


(km s-i) (K km s-i) 


12 CO J=3-2 


37.0 


67.0 92.4 ±1.9 


13 CO J=3-2 


44.5 


60.2 64.2 ±1.8 


HCO+ J=4-3 


45.5 


61.6 15.5 ±0.7 


CS J=7-6 


47.7 


60.7 16.6 ±0.8 




20 



-20 
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Figure 6. HCO + J=4— 3 emission integrated between 45 and 62 
km s — . The beam size is 22". 




10 
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Figure 7. CS J=7-6 emission integrated between 45 and 62 km 
s~l. The beam size is 22". 



peaks at a velocity close to that of the HCO + dip, we con- 
clude that this emission indeed appears self-absorbed. This 
kind of spectral feature might be revealing the existence of a 
density gradient in the clump (Hiramatsu et al. 2007), con- 
sistent with the presence of an embedded central source with 
outflowing activity. It is known that such molecular species 
enhances in molecular outflows (Rawlings et al. 2004). In 
effect, a strong enhancement of the HCO + abundance is ex- 
pected to occur in the boundary layer between the outflow 
jet and the surrounding molecular core. This would be due 
to the liberation and photoprocessing by the shock of the 
molecular material stored in the icy mantles of the dust. 
Figure 6 displays the HCO + J=4-3 emission integrated be- 
tween 45 and 62 km s _1 , showing a HCO + molecular clump 
peaking at the position of EGOg35. 

It is worth noting that the central HCO + and CS spec- 
tra, which survey the densest region where EGOg35 is em- 
bedded, have the closer component (i.e. the component with 
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lower velocity, or blue component) stronger than the farthest 
one (red component). This signature suggests the presence 
of infalling gas, which is consistent with the presence of a 
YSO accreting material. As Zhou et al. (1993) explain, this 
effect reflects the fact that the excitation of molecules is 
higher in the cloud center and, if the cloud has infalling mo- 
tion, then the observer is looking at the hotter side of the 
blue hemisphere and the cooler side of the red hemisphere. 
Therefore, the blue emission should always be as strong as, 
or stronger than, the red emission. 

Finally, it is important to note that the CS J=7-6 emis- 
sion maps the dense envelope where the massive YSO is 
evolving. The detection of this line implies the presence of a 
gaseous envelope with temperatures and densities above 40 
K and 6x 10 6 cm~ 3 , respectively (e.g. Takakuwa et al. 2007). 
Figure 7 shows the CS J=7-6 emission integrated between 
45 and 62 km s" 1 . 

3.1.1 Column densities and abundances 

To estimate the molecular column densities and hence the 
abundances in the region we assume as a first order approxi- 
mation local thermodynamic equilibrium (LTE) and a beam 
filling factor of 1. From the peak temperature ratio between 
the CO isotopes ( 12 T m6 / 13 T m6 ; taken from the CO mam 
components at v ~ 54 km s _1 ), it is possible to estimate the 
optical depths from (e.g. Scoville et al. 1986): 



13T 

1 Tl 



1 - exp(-T 12 ) 
1 - exp(-T 12 /X) 



(1) 



where t\i is the optical depth of the 12 CO gas and X = 
[ 12 CO]/[ 13 CO] is the isotope abundance ratio. Assuming 
R = 8 and using [ 12 CO]/[ 13 CO] = (6.21 ± 1.00)D G c + 
(18.71 ± 7.37) (Milam et al. 2005) where D G c = 5.5 kpc is 
the distance between the source and the Galactic Center, 
we obtain [ 12 CO]/[ 13 CO] = 53 ± 12. Then the 12 CO J=3- 
2 optical depth is t 12 ~ 300, while the 13 CO J=3-2 optical 
depth is T13 ~ 6, revealing that both lines are optically thick. 
Thus, we calculate the excitation temperature from 



T ex (3^2) 



16.95K 



(2) 



ln[l + 16.59K/(T max ( 12 CO) + 0.036K)] 

obtaining T ex ~ 20 K. Then, we derive de 13 CO and 12 CO 
column densities from (see e.g. Buckle et al. 2010): 



N( 13 CO) = 8.28 x 10 13 Tex + °- 88 



1 - exp(^L) 



15.87 s / r i3 dv ( 3 ) 



and 



N( 12 CO) = 7.96 x 10 13 ^ + °jg / r 12 dv (4) 

1 exp( T ^ ) J 

where, taking into account that r > 1 in both cases, we use 
the approximation: 



/ 



r dv = 
with 



1 



1-C-T / 



J(T)= , hv 



J(T ex ) - J(Tbg) 
hv/k 



T mb dv 



exp(^) - 1 



(5) 



(6) 



We obtain N( 13 CO) ~ 7x 10 16 cnT 2 and N( 12 CO) - 1 x 10 1 



Additionally, we derive the CS J=7-6 and HCO+ J=4- 
3 column densities from: 



N(CS) = 2.54 x 10 11 T - + °- 56 



1 — exp(- 



rdv 



(7) 



and 



N(HCO+) 



5.85 x 10 10 e 



25/l T ex + 0.71 



exp(- 



I 



rdv (8) 



Since the CS J=7-6 line is optically thick (e.g. Giannini 
et al. 2005), we use again the approximation presented in 
equation (5) assuming res = 1- On the other hand, by as- 
suming that the HCO + J=4-3 is optically thin, we use the 
approximation: 



/ 



t dv — 



J(T ex ) - J(Tbg) 



dv 



(9) 



We adopt as excitation temperatures 66 K and 43 K for the 
CS and HCO + , respectively, which correspond to the equiv- 
alent temperature of each molecular transition. We therefore 
obtain N(CS) ~8x 10 13 cm~ 2 and N(HCO+) ~ 9.5 x 10 12 
cm~ 2 . 

Assuming the abundance ratio [H 2 ]/[ 13 CO] = 77 x 10 4 
(Wilson & Rood 1994), from the N( 13 CO) we can estimate 
the H 2 column density in N(H 2 ) ~ 5 x 10 22 cm~ 2 . This 
value is in agreement with the H 2 column densities reported 
for high-mass protostar candidates associated with methanol 
masers (Codella et al. 2004; Szymczak et al. 2007), as in 
our case. Thus, using this H 2 column density we derive the 
abundance ratios for the HCO+ and CS: X(HCO+) ~2x 
10~ 10 and X(CS) ~ 1.6 x 10~ 9 . The X(CS) value is within 
the very wide range (5 x 10~ 10 — 2 x 10~ 7 ) measured in 
the outflows of protostars (e.g. Bottinelli & Williams 2004; 
J0rgensen et al. 2004; Giannini et al. 2005). 



3.1.2 High velocity material and outflows 

As discussed in Section 3.1, the 12 CO J=3-2 spectrum ob- 
tained towards the center of EGOg35 shows several com- 
ponents, some of them suggesting outflowing activity from 
EGOg35. The presence of outflows or high velocity material 
moving along the line of sight can be proven by compar- 
ing the 12 CO emission with the higher density tracer CS 
J=7-6. Figure 8 presents the 12 CO and CS spectra towards 
the center of the region with vertical lines remarking the 
ranges where 12 CO emission is detected at higher and lower 
velocities with respect to the CS. This figure confirms the 
presence of 12 CO spectral wings, which should be due to 
a red outflow going from 59.6 to 66.5 km s _1 and a blue 
one, extending from 37.0 to 48.9 km s~ x . In the case of the 
blue wing, we consider that the weak 12 CO component cen- 
tered at ~41 km s _1 is part of the source outflows because 
its velocity is compatible with that of the 6.7 GHz methanol 
maser detected by Cyganowski et al. (2009), confirming that 
this 12 CO component is due to the gas expansion caused by 
EGOg35. 

To analyze the velocity and spatial distribution of the 
12 CO J=3-2 emission, in Fig. 9 we present integrated veloc- 
ity channel maps every ~1.1 km s _1 . The spatial distribu- 
tions of the blue and red wings are shown in channels going 
from ~ 36 km s _1 to 50 km s _1 , and ~ 57 km s _1 to 64 
km s _1 , respectively. Figure 10 displays the 4.5 /xm emission 
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Figure 8. 12 CO J=3-2 (solid line) and CS J=7-6 (dashed line) 
towards the center of EGOg35. The vertical lines remark the 
ranges where 12 CO emission is detected at higher and lower ve- 
locities with respect to the CS. 



with contours of the 12 CO J=3-2 line integrated from 37 to 
49 km s -1 and 58 to 66 km s _1 (blue thick and red thin 
contours, respectively). 

In what follows, from the 12 CO spectral wings we es- 
timate the mass and energy to study the dynamics in- 
volved in the EGOg35 outflowing activity. Using M = 
fj.mHd 2 QX(CO)~ 1 N(CO) 1 we obtain the mass for the red 
and blue molecular outflows. N(CO) is the 12 CO column 
density, d the distance, mu the hydrogen atom mass, we 
adopt a mean molecular weight per H2 molecule of fj, — 2.72 
to include helium, Q. is the area of the 12 CO red and blue 
clumps shown in Figs. 9 and 10, and X(CO) = 10 -4 is 
the 12 CO relative abundance (Frerking et al. 1982). Inte- 
grating the 12 CO emission from 58 to 66 km s _1 we ob- 
tain a N(CO) red ~ 2.1 x 10 17 cm -2 and M red ~ 5 M , 
while integrating from 37 to 49 km s _1 we obtain N(CO)w ue 
~ 1.1 x 10 18 cm" 2 and M blw ~ 24 M Q . We calculate the 
momentum and energy of the red and blue components us- 
ing: 



P = MV 



E k 



1 9 
-MV 2 
2 



(10) 
(11) 



where V is a characteristic velocity estimated as the dif- 
ference between the maximum velocity of detectable 12 CO 
emission in the red and blue wings respectively, and the 
molecular ambient velocity (~ 51.5 km s _1 ), being V re d ~ 
14.5 km s~ and Vuue ~ 14.5 km s . Thus, we obtain 
Pred = 72 M Q km s" 1 and E r k ed = 510 M Q [km s" 1 ] 2 
(E r k ed ~ 1 x 10 46 ergs), and P Uue = 350 M km s" 1 and 
E b k lue = 2500 M Q [km s" 1 ] 2 {E b k lue ~ 5 x 10 46 ergs). The 
obtained outflows parameters are summarized in Table 3. 

The derived mass and energy for the outflows discovered 
towards EGOg35 are similar to those of massive and ener- 
getic molecular outflows driven by high-mass YSOs (Beuther 
et al. 2002; Wu et al. 2004). We do not provide estimates 
of outflow dynamical timescale and mass rate because the 
analysed high velocity material is aligned along the line of 
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Figure 9. Integrated velocity channel maps of the 12 CO J=3— 2 
emission every ~ 1.1 km s _1 . The grayscale is displayed at the 
top of each figure and is in K km B , the contour levels are 0.05, 
0.2, 0.4, and 0.7 K km s" 1 for the first panel, 0.5, 1, 3, 4, 4.5, 10, 
11, 12.5 K km s _1 for the second panel, and 0.5, 1, 2.5, 5, 7.5, 
10, 12.5, 15 K km s _1 , for the third one. The beam size is 22". 
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Figure 10. In grays is displayed the 4.5 /im emission from 
EGOg35. The blue thick and the red thin contours show the 12 CO 
J=3— 2 emission integrated from 37 to 49 km s — 1 and 58 to 66 km 
s — x , respectively. The levels are 6.5, 9.5, and 12.5 K km s — 1 for 
the blue emission and 4, 10, and 16 K km s^ 1 for the red one. 
The dashed square is the surveyed area. 



Table 3. Outflow parameters. 



Shift N(CO) MP E k 

(xlO 18 cm" 2 ) (M Q ) (M Q km s" 1 ) (Mourns" 1 ] 2 ) 



Red 0.21 5 72 510 

Blue 1.10 24 350 2500 



sight, and therefore it is not possible to estimate the lengths 
of the outflow lobes. In fact, with the moderate angular res- 
olution of the present observations, it is possible that we 
miss the emission of the outflows along the plane of the sky. 
Molecular observations with higher angular resolution are 
needed to spatially resolve the outflows. 
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Figure 11. Spectral energy distributions of the EGO 
G35.03+0.35. The black line correspond to the best fitting model 
and the gray lines to the following 19 best models. The dashed 
line is the photosphere model of the central source in the best 
fitting model. 



3.2 Spectral energy distribution 

Petriella et al. (2010) performed a spectral energy distribu- 
tion (SED) of EGOg35 concluding that the source is indeed 
a massive YSO. Taking into account that at present it is 
available new mid-IR fluxes obtained from WISE Prelimi- 
nary Release Source Catalog 1 , in this section we perform a 
new SED of the source to study with more accuracy its phys- 
ical parameters. The advantage of this new SED is that the 
WISE data present fluxes at 3.4, 4.6, 12, and 22 /im bands, 
being of importance the last two bands, mainly the 12 /im, 
because the SED of a YSO usually presents the separation 
between the contributions from the disk and envelope fluxes 
around this wavelength. Thus, we fit the SED using the tool 
developed by Robitaille et al. (2007) 2 . We use the fluxes 
in the four Spitzer-IRAC bands obtained from Cyganowski 
et al. (2008) together with the WISE fluxes, and the SCUBA 
bands at 450 and 850 fim (source G35.02+0.35, Hill et al. 
2006). The fluxes at 22 /mi from WISE and SCUBA bands 
were taken as upper limits as they have lower angular reso- 
lution and they can include contributions from other sources 
around EGOg35. As done in Petriella et al. (2010) we assume 
an interstellar absorption between 15 and 35 magnitudes and 
a distance range between 3 and 4 kpc. In Fig. 11 we show 
the SEDs of the 20 best fitting models. The solid black line 
represents the best fitting model from which we obtain the 
following parameters: M* ~ 23 Mq, M env ~ 1280 M@, and 
M env ~ 0.03 M-a/yr. We do not report the disk mass be- 
cause is not well constrained in the SED fitting. As Robitaille 
et al. (2007) point out, in the early stages of evolution, when 
the disk is deeply embedded inside the infalling envelope, the 
relative contributions of the disk and envelope to the SED 
are difficult to disentangle. 

Following the Robitaille et al. (2006) evolutionary clas- 
sification, the best fitting model (and also the following 80 
fitting models) correspond to a stage I YSO, i.e., a proto- 
star with large accretion envelope. In Fig. 12 we present the 
histograms with the distribution of the constrained physi- 
cal parameters. The hashed columns represent the 20 best 
fitting models and the gray columns correspond to all the 
models of the grid. We remark that the 20 best fitting mod- 
els correspond to massive central sources (between 8 and 
20 Mq) surrounded by massive envelopes with large accre- 
tion rates. Regarding the age of the source, there is a high 
dispersion in the results (ages range from 10 3 to 10 5 years) 
but the presence of a massive envelope is a strong evidence 
that the source is at the earliest evolutionary stage. 



3.3 Radio continuum analysis 

We investigate the radio continuum emission at 8.4 GHz of 
EGOg35. To produce the image we used archival data from 
observations performed with the VLA telescope operating 
in its B configuration on May 7 and 14, 2009 (project code 
AC948). Data processing was carried out using MIRIAD 



1 The Wide-field Infrared Survey Explorer (WISE) is a joint 
project of the University of California, Los Angeles, and the Jet 
Propulsion Laboratory/California Institute of Technology, funded 
by the National Aeronautics and Space Administration. 

2 http:/ /caravan. astro. wisc.edu/protostars/ 
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Stellar age (yr) 

Figure 12. Distribution of some physical parameters (central 
source mass, envelope mass, and age) of the 20 best fitting models 
(hashed columns) together with the distribution of all the models 
(gray columns). 



to 

<N 



<N 



O 

<N 




01.4 



18:54:00.7 



53:59.9 



Figure 13. Radio continuum emission at 8.4 GHz displayed in 
contours of 0.3, 1, and 10 mjy beam -1 over the TRAC-Spitzer 
4.5 /^m emission. The radio continuum has a beam of 1''4 X l'.'O 
while the angular resolution of the mid-IR emission is about 1'.'5. 



software package. The final image has an angular resolu- 
tion of l'.'4 x l'.'O and an rms noise of 0.1 mjy beam - . 
Figure 13 shows, in contours, the radio continuum emission 
of EGOg35 over its mid-IR emission at 4.5 /im. As can be 
seen from this figure, there are three radio sources: source 
A which coincides with the IR southwestern lobe of the 
EGOg35, source B that appears in between the IR lobes, 
and source C, marginally detected, which is projected onto 
the northeastern lobe of the EGO. The radio source A, the 
brightest and the only one resolved with this data set, is cen- 
tered at 18 h 54 m 00.49 s , +02°01'18.20" (J2000), and has an 
integrated flux density of 13.5 mjy. The source B is cen- 
tered at 18 h 54 m 00.66 s , +02°01'19.40" (J2000) and has a 
peak intensity of 1 mjy beam -1 , and source C is centered 
at 18 h 54 m 00.76 s , +02°01'22.62" (J2000) and its peak has 
a flux density of 0.3 mjy beam -1 . Source A and B were 
also detected at 44 GHz with a similar angular resolution 
by Cyganowski et al. (2009). The authors reported an inte- 
grated flux density of 12.7 mjy for the brightest radio source 
and a peak flux density of 3.6 mjy beam -1 for source B. We 
calculate the spectral index a (S oc u a ) between both fre- 
quencies, being a ~ —0.05 for source A, and a ~ +0.77 for 
source B. The radio spectral index of the radio source A is 
consistent with emission from an UCHII region, while the 
value obtained for source B suggests that its radio emission 
could be due either to an HCHII region (Kurtz 2005) or 
to a constant-velocity ionized wind source (Panagia & Felli 
1975). 

From this radio continuum analysis we can suggest the 
presence of several young stellar objects in the region. The 
fact that these radio sources overlap the 4.5 fim emission re- 
inforces the hypothesis that the infrared emission originates 
in YSO shocks. 



4 SUMMARY 

The extended green object EGO G35.03+0.35 (EGOg35, 
in this work), a massive YSO, is embedded in a molecular 
clump located at the border of an HII region. In this work we 
investigated the surrounding molecular gas through several 
molecular species using the Atacama Submillimeter Tele- 
scope Experiment (ASTE). We observed and analysed the 
12 CO and 13 CO J=3-2, HCO+ J=4-3 and CS J=7-6 tran- 
sitions, which are useful to trace outflows and dense gas. To 
complement our analysis we used IR and radio continuum 
data from public database. In what follows we summarize 
the main results of our work: 

(1) Most of the molecular spectra observed towards the 
surveyed region are far of having a simple Gaussian shape, 
presenting asymmetries, absorption dips, and spectral wings 
or shoulders, characteristics that strongly suggest the exis- 
tence of kinematical perturbations in the gas originated in 
the object EGOg35. 

(2) The 12 CO J=3-2 line towards EGOg35 shows a 
double peak structure like with a main component centered 
at ~54 km s , a less intense component centered at ~48 
km s , and an absorption dip at ~51.5 km s . The line 
also presents spectral wings due to the ouflowing activity of 
EGOg35. We obtained a total mass and kinetic energy for 
the YSO outflows of 30 M and 3000 M [km s -1 ] 2 (6 x 10 46 
ergs), respectively, similar to those found towards other mas- 
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sive and energetic molecular outflows driven by high-mass 
YSOs. 

(3) We discovered a HCO + clump towards EGOg35, 
supporting the presence of outflows in the region because it 
is known that such molecular species enhances in the bound- 
ary layer between the outflow jet and the surrounding molec- 
ular core. The central HCO + spectra present an absorption 
dip at v ~ 54 km s _1 due to self-absorbed gas, which is 
probably revealing the existence of a density gradient in the 
clump. 

(4) We discovered a CS clump towards EGOg35. The 
CS J=7-6 emission maps the dense envelope where the mas- 
sive YSO is evolving, and its detection implies the presence 
of molecular gas with temperatures and densities above 40 
K and 6 x 10 6 cm -3 , respectively. 

(5) The central HCO + and CS spectra, which survey 
the densest region where EGOg35 is embedded, show two 
velocity components with the closer one (blue component) 
stronger than the farthest one (red component). This sig- 
nature suggests that these lines are mapping infalling gas, 
which is consistent with the presence of a YSO accreting 
material. 

(6) From the spectral energy distribution (SED) study 
we confirm that EGOg35 is a massive YSO at the earliest 
evolutionary stage (i.e. a class I YSO). 

(7) Analysing radio continuum emission towards 
EGOg35 we conclude that there is evidence of the presence 
of an UCHII region and another source that could be ei- 
ther an HCHII region or a constant-velocity ionized wind 
source. The radio continuum analysis suggests the presence 
of several possible young stellar objects in the region. Fu- 
ture multiwavelength observations with very high angular 
resolution will be of a great importance to go deeper in the 
study of this region. 
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